Angular distributions in hard exclusive production of pion pairs 
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Using the leading order amplitudes of hard exclusive electroproduction of pion pairs we have 
analyzed the angular distribution of the two produced particles. At leading twist a pion pair can be 
produced only in an isovector or an isoscalar state. We show that certain components of the angular 
distribution only get contributions from the interference of the 7=1 and the (much smaller) 1 = 
amplitude. Therefore our predictions prove to be a good probe of isospin zero pion pair production. 
I ' We predict effects of a measurable size that could be observed at experiments like HERMES. We 

also discuss how hard exclusive pion pair production can provide us with new information on the 
f^*) ■ effective chiral Lagrangian. 

<n : 
x> : 

r*"- ■ I. INTRODUCTION 

The factorization theorem jjj states that the amplitude of the exclusive production process 

oo ! 

H + T -» F + T> (1) 

' at large invariant collision energy ^fs — > oo, large virtuality of the (longitudinally polarized) photon Q 2 — ► oo, fixed 
7-H Q 2 / s (Bjorken limit), and with — t, m 2 ^, mj./, mr F <SC Q 2 can be written in the form 
O; 

2_. / dzdxifE T (xx, xi — xbj, t) Hij ( ,Q 2 ,z) &J(z) + power-suppressed corrections , (2) 

q-: 

where is a hard part computable in pQCD as a series in a s , <&J is the distribution amplitude of the hadronic state 
F, and fTL, is a T — ► T 1 skewed parton distribution (for a review see ||). Skewed parton distributions (SPD's) 
' are related to matrix elements of bilocal operators between states of different momentum. They are generalizations 
^ , of the usual parton distributions which parameterize the diagonal matrix elements of the corresponding operators. 
Therefore hard exclusive electroproduction opens a new way to study the partonic structure of hadrons. Naturally, 
for the experimentally accessible range of Q 2 the size of the 1/Q suppressed power corrections can be rather large, 
the quantitative estimates of higher twist corrections are therefore in the focus of recent investigation & . Still, the 
formalism has been used to investigate exclusive production of single light mesons |p],[7|-|l0|] . For the starting point of 
the present work it is important to note that the factorization theorem is not limited to the case where the produced 
hadronic state F in reaction (Q) is a one particle state provided that its invariant mass is small compared to Q [[ll| . 

In the present paper we present a detailed analysis of hard exclusive electroproduction of two pions to leading 
order in the strong coupling constant a s . In that case $^ in the expression (|J) is a two-pion distribution amplitude 
(27rDA). At leading twist the pion pair can be produced in a state with isospin one or zero. For 1 = 1 the process is 
dominated by the p peak. (So in a certain sense our approach can be considered as an alternative description of p 
production, which has been studied earlier using a distribution amplitude for thep meson @,| 0). The description 
of pion pair production at small xej in terms of 27rDA's was considered in Ref. 113] , where it was demonstrated that 
such a description allows us to determine the details of the partonic structure of the pion and the p meson from data 
on the di-pion mass distributions in hard diffraction alone. 

Here we are mainly interested in the production of / = pion pairs. It seems to be difficult to investigate isoscalar 
pair production by the measurement of total cross sections because identifying 7r°7r° production events would require 
detecting 4 correlated photons whereas 7r + 7r~ production is dominated by the isovector contribution due to the p 
resonance (see our analysis in [p"3|). However, the angular distribution of the produced pion pairs contains components 
that depend only on the interference of the 1 = and the 1=1 channel Jl3],[l4| . Using specific models for the involved 
SPD's and 27rDA's we predict effects of a measurable size and suggest the most favorable kinematic range for such 
measurements. 
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The interesting feature of hard exclusive production of pion pairs in an isoscalar state is that the pions are produced 
not only by a collinear qq pair but also (and in the same order in 1/Q and a s ) by two collinear gluons. This observation 
opens a possibility to access the gluon content of the isoscalar irir states. In the analysis of the present work we shall 
assume that the elusive isoscalar tttt resonance Jo (400 — 1200) has no anomalously large gluon contents as it was 
suggested e.g. in [ ^5|Jl^ ]. The use of hard exclusive pion pair production for "gluonometry" of the low-lying isoscalar 
resonances will be considered elsewhere. 



II. AMPLITUDES AND INTENSITY DENSITIES 



In this section we calculate the amplitude of hard two-pion electroproduction to leading order in 1/Q 2 and a s in 
terms of SPD's and 27rDA's and show how they are connected to the so called intensity densities which we define as 



(Pi (cos 0))** 



j\d cos Pi (cos 9] 
/^dcosfl d<7 " 



d cos 9 



(3) 
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k + k = q' 



FIG. 1. Definition of the scattering angle 8. 




Here 8 is defined as the scattering angle of the ir + or one ir° (for it + it or ir°ir° production respectively) in the center 



of mass system of the two pions with respect to their total momentum e.g. in the laboratory system (see Fig. |l|)[|]. 
We consider the process 

1 Uq) + B 1 (p)^7TTT(q') + B 2 (p + A), (4) 

where a linear polarized virtual photon with momentum q and a nucleon B± with momentum p produce a final state 
with a baryon B 2 (a nucleon or a heavier state) with momentum p' = p + A and a pion pair with total momentum 
q' . Note that the amplitude for the corresponding production process with a transversally polarized virtual photon 



is 1/Q power suppressed (with Q = —q ). Therefore the (sub-) process 
is the only leading contribution to the electroproduction reaction e~(l) + 
Using two light cone vectors n and h normalized such that 



n-(p+j/) = 2, 
the relevant momenta can be expressed as 



n ■ n = 1 



jj) with a longitudinally polarized photon 
Bi(p) -> 7T7r(g') + e -(l-q) + B 2 (p + A). 



(5) 



p" = (l + 0#i" + (1-0^-^ 

1 



2 " '2 x 



o 2 

q* 1 = -2£n^ + -^n M 
4£ 



(6) 



with 



1 In principle one can analogously define a more general intensity density {Yi m (6, 4>)) 7Tn . However, our analysis shows that at 
leading twist the production rate is independent of the azimuthal angle cj>. Therefore non-vanishing results are obtained only 
for the intensity densities (Yio(0, 0))*^, which are up to a constant normalization factor equal to (Pi(cos 6)) 717 ' defined in (H). 
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1/ 2 

e = ^l-i(?' 2 -*-4e 2 rn 2 ))+0(i i ). (7) 



i / o 2 ^ 

= - (m R , + m Ba - - 



Aj_ is defined to be that spacelike component of the momentum transfer A which is transverse to the light cone vectors 
n and n, t is the momentum transfer squared t = A 2 , and £ is the skewedness parameter describing the longitudinal 
component of the momentum transfer 

e = ~(«-A), (8) 



which in the Bjorken limit Q 2 — > oo can be expressed in terms of the Bjorken variable xbj = ^ 

2 - I B j ' 



(9) 



S| B2 5| fi. 5| fi. 

FIG. 2. Examples of the three classes of leading order diagrams for hard exclusive pion pair production. 

The leading order amplitude for the reaction (Q) corresponds to diagrams of the type shown in Fig. || and has the 
form (see also @@) 



(B 2 (p'),n a n b ( q >)\r- m - -SlIB^p)) 
-{eA*a.)± f dz 



Cf 



Y,F ff ,(T,t,t)<P¥ f (z,(,m^) 



/,/' 



z(r + C-ie) (l-z)(r-£ + ie) 
1 1 



t + £ - ie r — £ + ie 
1 1 



2N C *—~ I 13 z(l — z) [T + £ — le t — £ + ie 



(10) 



where e/ is the electric charge of a quark of flavor / in units of the proton charge e and is the longitudinal photon 
polarization vector 



, Q. 



W H — =rr 



Q 4£ 



(11) 



(=> = i(g 3 ,0,0,9°) in a frame with n" oc (1,0,0,1) and oc (1,0,0,-1)). The functions F Sf ,{ T ,^t) and 
_Fg(t, £,£) are skewed parton distributions for quarks or gluons respectively defined as 



*fe(r,e,t) = ± / ge^W(B 2 (p')|T[^(-^)^(f )]|B 1 (p)> 
The gluon color index A is understood to be summed over. The two-pion distribution amplitudes 



'An' 



(12) 
(13) 
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*&(z,C,m™) := / ^e- iA ^«'^( 7 r o 7r b |T[v; / (An)^/'(0)]|0) (14) 



and 



^(^C,m») := ^e- iA ^'«)^^(7rV b |T[i^«(Afx)F^(0)]|0) (15) 

describe how a pion pair is produced by two quarks or gluons respectively. The variable z is the fraction of the 
longitudinal (along the vector n) meson pair momentum q' = q — A carried by one of the two partons. The variable £ 
characterizes the distribution of the longitudinal component of q' between the two pions and is given in terms of the 
momentum k a of the pion 7r a and the total momentum of the pion pair q' by 

, k a ■ h 

C=^^- 16 
q ■ n 

Let us note that £ is related to the angle 6 defined above by 



/ 4tti 2 

/3cos0 = 2(-l with 0~\l ^ . (17) 



In the definitions ( |l2| ) - (15) a gauge link of the form V[exp(ig j* 1 A^dz^)] with = ^2^t A A A is implied to be 
inserted between the two operators at different space-time coordinates. Also we do not write out explicitly the scale 
dependence of the SPD's and 27rDA's. 

In the following we will restrict ourselves to the case where B\ and B2 both are proton states. Then the pion pair 
can be either 7r + 7r~ or 7r°7r°. The 2ttDA's @ and © can be isospin decomposed and expressed in terms of only 
two independent quark 27rDA's corresponding to isospin 1 = and 1=1 pair production respectively and one single 
gluon 2-7rDA. For 7r + 7r~ production we have 

$J, + />,C,to™) = 8 f ' f * I=0 {z,C,m mr )+>i{' f * I=1 {z,C,m mr ) (18) 

$£ +7r ~ (z, C, m ffjr ) = $ G (z, C, m wv ) (19) 

and for 7r°7r° production 

l},°/feC,m II )=A I = ,) (2,C ) m II ) (20) 

^(«,(,m„) = ^ G {zX,m^). (21) 

Taking into account the C-invariance of the underlying theory one can easily derive the following symmetry properties 
for the three 27rDA's: 

$ 7=0 (z, £ m™) = -$ 7=0 (1 - z, C, m ffjr ) = $ /=0 (z, 1 - C, to^) (22) 
$ 7=1 (z, C, m^) = $ /=1 (l - z, C, to™) = 1 - C, to™) (23) 

$ G (z, C, to™) = $ G (1 - z, C, m™) = $ G (z, 1 - C, m™) (24) 

Using the decompositions (18) - ( f2l| ) and the symmetry relations (p2|), ( p3| ) the transition amplitudes ( fTo| ) for 7r + 7r~ 
and 7r 7r production and an elastically scattered proton target can be written in the form 

T n+ ^ = T I=Q + T I=1 (25) 
r V = T I=0 (26) 



with 



-.7=1 _ /„,i 1 5,7 = i/ nT + , r+ , 3 



- -( e4 ™^i^* l 27 - +I * + ^ • (28) 



Here $ 7_0 and are functions of m™ and cos# defined as 
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7=0 



dz- 



z(l-z) 



7=1 



= / dz 



2(1"*) 



and , /g are the following integrals over SPD's depending on xej and t: 

1 1 



^ = f drF ff (T,Z,t) 
Ig = J drF G (r,e,i) 



t + £ — ie r - £ + ie 
1 1 



t + £ — ie r - £ + ie 



(29) 
(30) 

(31) 
(32) 



In the next section we will show that restricting the ^-dependence of the 27rDA's to their asymptotic shapes but 
keeping the momentum fraction carried by quarks in the pion (gluons have the momentum fraction Mp = l — M®) 
as a free parameter these functions take the form 



$ J = u (z,C,m™) 



120M? , 

2 -z(l-z)(2z-l) 



3C — 1 1 

— ^ — fa(m ww )P (cos6) - — / 2 (m W7r )P 2 (cos( 



$ I=1 (z, C, m^) = 6z(l - «)(2C - l)i^(m w ) 
$ G (z, C, m X7r ) = -60M 2 G z 2 (l - zf 



3C — (3 2 (3 2 
^2 — /o("i 7r7r )P (cos6») - — / 2 (m W7r )P 2 (cos( 



(33) 
(34) 
(35) 



where /o and / 2 are Omnes functions, F^im^) is the electromagnetic form factor of the pion, Nf the number of light 
quark flavors, and C = 1 + fern 2 with b = — 1.7GcV~ 2 is an integration constant estimated in the instanton model of 
the QCD vacuum [jl7|. Using these expressions the integrals ( p9| ) and ([30j) can be evaluated: 



40M, 



20M 2 G 



3C — /3 2 /3 2 

^ — fQ{m^)P Q {cos6) - — / 2 (m Tr7r )P 2 (cos6») 



i.'- 1 = 6/3P 7r (m W7r )Pi(cos( 



(36) 
(37) 



The differential cross section of the electroproduction process is given in terms of its T-matrix element T" n by 

da*""*" m^P 



dxBjdydt dm-xTT d cos 8(4-7r) 5 p • Z 



E 



(38) 



where I is the initial lepton momentum, y is defined by y = p-q/p-l and corresponds to the energy loss of the scattered 
lepton in the proton rest frame, and the sum is understood to run over the spin polarizations of the scattered nucleon 
and lepton. The T-matrix element of the electroproduction reaction e~(l) + B\ (p) — > Tm(q') + e~ (I — q) + P 2 (p + A) 
is related to the amplitude of the sub-process (0) calculated above by 



= u(l - q, S ')^u(l,s)e^T^ . 



For the spin sum one gets 



El* 



7T7T I 2 



2e\l-y) 



E' T 



7T7T I 2 



S's' 



p-lxBjy 3 ^7 

Using ( |38| ) and fllcf ) the intensity densities defined in (||) are obtained to be 



(39) 



(40) 



(P(cos#))™ = 



J dx Bj dydtdm^ m„ T /3{l - y)y- 3 x B ^ j\ d cos 0P ; (cos 0) £ s , l T ™l 2 



J dxBjdydtdm^ 771,^/3(1 - y)y 3 x Bj 1 J_ x d costs' \ T ™ 



(41) 



The cos 0-integrals can be solved analytically, they only involve integrals over products of 3 Legendre polynomials as 
can be seen comparing the previous equation with (p|) - (pi), (|36|), and (|37|). The full expressions are rather lengthy 
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and are shown in the appendix. Non-vanishing results are obtained for 7r + 7r~-production for I = 1,2,3,4. Note that 
for I = 1 and I = 3 only the interference term (^i=oyrpi=i _|_ cc contributes to the nominator of Eq. (f4l|). These 
two intensity densities are highly sensitive on the production amplitude for isoscalar pion pairs, a vanishing 1 = 
amplitude would imply that (Pi(cos0)) 7r+7r and (P3(cos^)) 7r+Tr are zero. For ^"^"-production we can predict that 
only the intensity densities (P2(cos8)) 7r * and (P^cos 6)}^ * do not vanish. 

We have mentioned above that our analysis is valid only up to 1/Q power corrections, which are not necessarily 
negligible at realistic scales in experiments like HERMES. Therefore it is desirable to define quantities for which at 
least some power suppressed contributions are absent. Indeed we can exclude the basically uncontrolled contribution 
of transversally polarized photons taking an appropriate linear combination of (Pi(cos0)) Tr7r and (P3(cos6')) 7r7r . The 
analysis of |l8| shows that the combination (Pi (cos 0))™ + \/7 JZ{P$ (cos 0))*"* projects out the component of vanishing 
helicity of the final two pion state (see also appendix B of |l9|). Assuming s-channel helicity conservation we can 
conclude that only longitudinal polarized photons contribute to this component. 



III. TWO-PION DISTRIBUTION AMPLITUDES 



The two-pion distribution amplitudes defined i n (|14| ) and ( |15|) describe the fragmentation of a pair of collinear 
partons (quarks or gluons) into the final pion pair ]20|]. Some properties of the 27rDA's were already presented in the 
previous section (Eqs. ( flq ) - (^J))- In this section we want to discuss further constraints for these functions that 
allow us to construct realistic models for them. 

Following |l^,^l| we decompose both quark and gluon 2-7rDA's in conformal and partial waves (see For the 

quark 27rDA's the decomposition reads 



OO 71+1 



/=0 



(z, C, m™) = 6z(l - z) E B I n f°(m^)C 3 J 2 (2z - l^/T^C - 1)) 



$ i - 1 (*,C,m» T ) = 6«(l-*) 



n=l 1=0 
odd even 

oo n+1 

£ ]T B I n f(m vn )C 3 J 2 (2z - l)Pi(p-\2( - 1)) , 

n=0 1=1 
even odd 



while for the gluon 27rDA we have 

$ G (z,C,m„) = 30z 2 (l 



OO 71+2 

z f E E AZ{m^)Cl'\2z- l)P(r 1 (2C- 1)), 

even even 



(42) 



(43) 



(44) 



with the Legendre polynomials Pi{x) and the Gegenbauer polynomials C%(x). The expansion of the z-dependence 
in Gegenbauer polynomials is chosen such that under evolution of the 27rDA's the coefficients B^f°, B^f 1 , and 
are renormalized multiplicatively with mixing only between B^f° and A^_ ± ; ^2j,^3|. With the choice of Legendre 
polynomials for the ^-expansion the above decomposition is directly related to the partial wave expansion of the 
resulting production amplitude for pion pairs. 

The 27rDA's are related to the distribution amplitudes of a single pion ip v (z) by soft pion theorems of the form 

$ /=0 (z, C = 1, = 0) = <P I=0 (z, C = 0, m n7r = 0) = (45) 
$ /=1 (z,C= l.m™ = 0) = -$ / = 1 (z,C = 0,m 7r , = 0) = <p{z) (46) 
$ G (z, C = 1, = 0) = $ G (z, C = 0, = 0) = . (47) 

Additional constrains are provided by the crossing relations between the quark and gluon 27rDA's and the correspond- 
ing (skewed and forward) parton distributions in the pion. For the derivation of the crossing relations see [ |l7fl (for 
quark 27rDA's) and [|T) (for the gluon 27rDA). Expressed in terms of the coefficients B T nl , and A G ; they take the form 



Cln(0) 



2 2n+ 1 

3 n+ 1 



d ^" _1 ^E(^)+^>)) 



2 2n+ 1 



^n-2,n(0) 



3 n + 1 j 

4 2n+ 1 



&xx n -\£ + (x)-t + (x)) 



5(n+l)(n + 2) 



dx : 



(48) 

(49) 
(50) 



G 



where ql(x), and g v (x) are the usual quark and gluon distributions in the pion. Using these relations one can easily 
derive the normalization of the 27rDA's $ /=0 and $ at m^n = 0: 

J dz(2z-l)$ / =°(z,C,m^ = 0) = -4- M 2 Q C(1-C) (51) 



N f 

/ dz$ G (z,C,m o»,fG 
Jo 



= 0) = -2M^((1 - C) . (52) 

Here and M G are the momentum fractions carried by quarks or gluons respectively in the pion. The first moment 
of $ /=1 can be determined by crossing even for arbitrary values of m W7r and is expressed in terms of the electromagnetic 
pion form factor F^^m^^) in the following way: 

dz<S> I=1 (z, C, m^) = (2C - l)F«{m^) . (53) 



This is due to the fact that crossing relates the moment on the left hand side of Eq. (|53| ) to the first moment of the 
corresponding pion SPD which is given by the pion form factor. 

Before studying the m^-dependence of $ /=0 and <I> G we want to discuss the asymptotic shape of the 27rDA's. In 
the asymptotic limit only the coefficients B^q°, Bj^°, B^f 1 , Aq Q , and Aq 2 do not vanish. Therefore the asymptotic 
expression for the isovector 27rDA has the form 



$ 7=1 (z, C, m™) = 6z(l - z)(2C - (m ffff ) , (54) 



while taking into account the soft pion theorems ([45]); fl47|) and the normalization conditions (|5l|), ( p2|) we get for 
$ /=0 and $ at asymptotically large Q 2 and vanishing mass m vv 

i on yr Q 

$ /=0 (z, £ m_ = 0) = £p-*(l - z)(2z - 1)C(1 - C) (55) 

$ G (z, C, m wff = 0) = -60M G z 2 (l - z) 2 C(l - ■ (56) 

The explicit form of the Q 2 -dependence in LO implies that the linear combinations Aq — Cf^B[q° and Aq 2 — Cf 

die out at large Q 2 due to the mixing between the gluon and the quark 2-7rDA. This can be used to fix the asymptotic 

values for M 2 Q and M§ to 

M Q ( ™,) = _Nf_ (5?) 
M2G(asymp ,) = _4g_^ (5g) 

However, in order to be more general we do not restrict ourselves to these values but keep M 2 as a free parameter. 

The m^-dependence of the 27rDA's can be analyzed by means of dispersion relations jlTf]. A similar analysis can 
be found e.g. in Refs. |^5|,|6), where the effect of final state interactions on the decay of a light meson into two pions is 
investigated. In the region to 2 ^ < 16m 2 the imaginary part of a 27rDA is related to the pion-pion scattering amplitude 
due to final state interactions by Watson's theorem p7j| . For <f> /=0 (In the analysis of <!> G the line of argument is 
exactly the same) this relation can be written in the form 

Irn^nrV™)} = sk^KjJe^^JE^K,)* = tan(5 i °(m 2 r7r ))Re{^f °(m™)} , (59) 



where Sf are the two pion phase shifts. From Eq. (g9|) the following iV-fold subtracted dispersion relation for the 
coefficients B^m^) can be derived: 



B**(rn 1-V m '" d * B'=<Vfrt + m - r d Jan(^( S ))Re{^ r °(^)} 

B nl M - —^wBm (0) + — ds sN{s _ m ^_ it) ■ (60) 



Solutions of such dispersion relations were found long ago by Omnes pq | and can be written in the forn^] 



2 We give here the simplest form of the solution with only one subtraction and assuming that B^f°( m mr) has no zeros in 
the relevant m,r„- interval. The presence of a zero in B^ Q (m mt ) at some value of m,^ can be interpreted as a signature for a 
glueball in the nn channel. This case shall be considered elsewhere. 
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^r°(o) 



/i( m 7T7r) = CXp 



M?«r) + 



-Re 



00 ds W 

2 s(s - to 2 w - ic) 



4 m 



(61) 



where //(m^) are the so called Omnes functions. By derivation, the solution ( |6l|) is valid only for m 2 ^ < 16m 2 , 
however, the deviations due to inelastic pion-pion scattering are expected to be small up to rather high energies. In 
terms of the Omnes functions the (quasi-) asymptotic 27rDA's <£> /=0 and $ G take the form 



r-n, , \ 120M 9 Q . 



2V 



/ 



3C — 1 (3 2 
— 77; — foim^PoicosO) - — / 2 (m TTr )P 2 (cos( 
12 



* (^C.nw) = -60M 2 G z 2 (l - z) 



3C — /3 2 /3 2 

^2 — /o(m 7r7r )P (cos6») - — / 2 (m 7r7r )P 2 (cos6') 



(62) 
(63) 



The constant C in Eqs. ( |62[ ) and fl63| ) plays the role of an integration constant in the Omnes solution of the corre- 
sponding dispersion relation. From the soft pion theorem it follows that C — 1 + 0(m 2 ). Using the instanton model 
for calculations of B n i(m^ T ) at low energies |l7|j2"i| one finds the constant C to be equal to: 

C= 1 + bml + 0(m*) with fr«-1.7GeV~ 2 . (64) 

In the section [v] we shall discuss the relation of the constant C and the near threshold behavior of f b(?ri 7r7r ) to the 
effective chiral Lagrangian. The expressions (Ejl), (62), and ( |63] ) are exactly those we have used for the 27rDA's in the 
analysis of the previous section. 



IV. MODELING THE SKEWED PARTON DISTRIBUTIONS 



The skewed parton distributions defined in Eqs. ( |l2| ) and ([l3]) can be decomposed into different Lorentz structures. 
Adopting the notation of || we can write 



in n^A" 

F ff (r, e, t) = H f (r, t) u(p', S')Mp, S) + E f (r, £, t) u{p' , S')-^ u(p, S) 



2niN 
\(j A v 

F G (r, t) = Hq(t, £, t) u(p\ S')M Pl S) + E g (t 7 t) u(p> , S 1 )^ 



u(p, S) . 



(65) 
(66) 



The SPD's Hf(r, £, t) and Hq{t, £, t) are highly constrained by their forward limit whereas only little is known about 
the functions Ef(r,Z,t) and Ec(t, £,t) (however, it is clear that their contribution vanishes when the skewedness £ 
goes to zero). 

For our analysis we model the functions Hf(r, £, t) and Hq{t, £, t) using Radyushkin's double distributions 0. The 
SPD's H(r,£,t) are related to the so called nonforward parton distributions !F^{X,t) by: 



e(C-r)jf^ r 



riJ G (r,^,<) = - 



O(r + 0^ ( 



,VTTt.'J- e «-^v 1 + e . 



(67) 
(68) 



with £ = The skewed parton distributions are obtained from the double distributions by the integral 



? < {X,t) = 



dyF(X~(y,y;t) 



(69) 



For the double distributions we adopt now the ansatz suggested by Radyushkin Q F(x,y;t — 0) = q(x)ir(x,y) 
(for quark distributions) and F(x,y;t = 0) = x G(x) n(x, y) (for gluon distributions), where q(x) and G(x) are 
the ordinary quark and gluon parton distributions in the nucleon and it{x, y) is a profile function chosen to be 
ir(x,y) = 6 y(l — x — y)(l — x)~ 3 for quarks and ir(x, y) = 30 y 2 (l — x — y) 2 (l — a;) -5 for gluons. The (forward) parton 
distributions are modeled using the MRS(A') parameterizations of the reference |3^] at Q 2 = 4GeV 2 (see footnoted). 



We take the SPD's at this fixed scale. For evolution effects see e.g. 1311 



8 



For the i-dependence of the SPD's we use a factorized ansatz ||,[T(| motivated by sum rules that relate their first 
moments to form factors p]. The quark SPD's are constrained by 



j\rH f (r^,t)=F( /p (t) 



(70) 



with the Pauli form factor Fi of the proton with respect to the flavor /. The corresponding ansatz for the symmetric 
part of Hf(r, £, t) is 



Hf{r^,t)+H f {-T^,t) 



F( /P {t) 
F{ /p (0) 



(71) 



We adopt a parameterization for the nucleon form factors taken from the appendix of |32| and take into account the 
relations 



F u/p = 2p p + p n ) pa/p = 2p n + pp ( 



d/p _ 



(72) 



which are valid if the contribution of strange (and heavier) quarks to the nucleon form factor is neglected. Taking the 
first moment of t£/g(t, £,i) we can analogously motivate the model 



F a (0) 



with the gluon form factor of the proton defined by 



1 



(p'\G^(0)Gi(0)\ P ) = (2p^p v - ^g^m 2 N )F e (t) . 

For the gluon form factor we can adopt the parameterization of []33f : 

F e {0) 



Fe(t) 



(73) 



(74) 



(75) 



with a — 3 and /i 2 = 2.6 GeV 2 . The same i-dependence we assume also for the antisymmetric part of Hf(r,£,t) 
entering the integral if (see Eq. (^l))). The combination Hf(r,£,t) — Hf(—T,£,t) is not constrained by the sum 
rule (f70|), however, we expect a behavior similar to the gluon SPD due to the fact that H/(t, £, t) — Hf(—r, £, t) and 
Hq(t, £, ,t) mix when they are evolved. 

In p4j it was shown that the SPD's obtained by double distributions are not complete and that an additional term 
has to be added, the so called D-term. In our analysis we take into account this result by adding to the Radyushkin 
model result for Hf(r,^,t) a term estimated from the chiral quark-soliton model. We write 



^(Radyushkin model) 



H f (T,Z,0) = H} 
taking for the function D(x) the following numerical estimate 



(T,e s o)+e(e-M)ii>(r/o 



D(x) = -4 (1 - x l ) Cl ,2 {x) + 0.3 Cl' A {x) + 0.1 C 5 J/z (x) 



3/2 , 



t3/2. 



(76) 



(77) 



extracted from a calculation of the singlet quark SPD in the chiral quark-soliton model |3q | . The same term with the 
opposite sign and multiplied with the corresponding form factor is taken as model for Ef(T,£,,t) in order to satisfy 
the sum rule S 



J dTTj2[Hf(T,Z,t) + E f (T,Z,t) 



= independent of £ . 



(78) 



(Note that the part of Hf(r,£,t) resulting from a double distribution automatically satisfies Eq. (|7q).) 
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V. PROBING THE EFFECTIVE CHIRAL LAGRANGIAN IN HARD EXCLUSIVE REACTIONS 



In this section we consider the case when the produced pion pairs have an invariant energy close to the threshold 
m-n-Tr = Am 2 . In this case the dependence of the intensity densities on m„ is related to the effective chiral Lagrangian 
(EChL) describing the interaction of soft pions with gravity. Therefore one can use data on hard exclusive pion pair 
production to probe this yet unknown part of the EChL. 

Due to the QCD factorization theorem (Q) one has a well defined separation of the short and large distance parts 
of the interaction. Since the large distance parts of the process are defined in terms of hadronic matrix elements of 
QCD quark and gluon operators which are independent of the hard scale (up to a logarithmic scale dependence which 
is controlled by well known evolution equations) one can apply the methods of effective chiral perturbation theory to 
describe properties of these matrix elements in a model independent way. 

In the analysis below we make the assumption that in hard electroproduction of pion pairs the two pions are 
produced dominantly by the operators of the lowest conformal spin (index n in Eqs. ((42|), (f43|), and (44)). This 
assumption is justified for large Q 2 since the contribution of operators with higher conformal spin logarithmically 
dies out with increasing Q 2 due to their larger anomalous dimensions. For the 1=1 channel the lowest conformal 
spin corresponds to the vector current operator and for / = to the energy momentum tensor. Therefore using this 
assumption the constant C in Eqs. (62) and ( |(33| ) and the near threshold behavior of the functions F n (m n7r ), ^(m^), 
and /^(w^), which enter the expressions for the various intensity densities (see the appendix for a collection of 
formulae), are fixed by the EChL. 

Let us write the relevant terms of the fourth order EChL describing the interactions of soft pions with electromag- 
netic fields and with gravity [{36|]37f l: 



£ (4) =-LL 9 Tr F»"dMd v U 1 + L n RTr dMd^ + L 12 R, V ^ d^Ud^U^ 



L 13 R Tr 



iU + U ] i 



(79) 



Here F^ v is the field strength of the photon field, and R are the Ricci tensor and the curvature scalar of an external 
gravitational field, and U = exp(m a \ a / f v ) is the non-linear pseudo-Goldstone field. The ellipsis stands for the terms 
of the EChL that are not relevant for us here. Now we can use the results of the calculations in Refs. |36|,[57| to express 
the constant C and the near threshold behavior of the functions /o, and f 2 in terms of the constants Li in the 
EChL ([79|). Let us first introduce the scale dependent constants L\ (/z) which appear after proper renormalization of 
ultra-violet divergences we have 



Ll(fi) = L i ~T l A, 

r 9 



1 1 _ 2 
-, Tn--, r 12 _o, ria--, 



a = 



rf-4 



(4tt) 2 



(80) 
(81) 

(82) 



for further details see Refs. ]36| , p7| . Using these definitions and the results of chiral perturbation theory |3^j38J we 
obtain for F- K {m^T T ) 



Ami 



where the standard loop integral is defined by 



1927T 2 



'-K K {m^) 



? / 2 
™ 1 1 + 2 In 



1927T 2 



J\ X In (l -,(!-,)!<) 



-fl+ln-£ , (83) 



(84) 



The result for the constant C and the near threshold behavior of f^m^^) and jiim^^) can be easily extracted from 
the results of Ref. Wh: 



C 



= 1 - 

4 



8(iL r n (fi) + L r M - 2L r 13 (fj,)) + -1- (5 In -4 ~ ^ ^ -4 + 2 In - ^ 

^ UAfV 12VfV LA\r-JJ 48^2 V m 2 3 m 2 m 2^ ^ 

(AW - /o(m™)) = 1 + 2^{4L r M + UM) - 16^| (Z^OO - L r M 



2m 



2 

7T7T 



16^ 2 /2 



—K n (mim) 



167T 2 /- 



f 2 



A8ir 2 f 2 



ji(jn^) = 1 - 2L\ 2 (fj,) '- 



f 2 



(85) 

(86) 
(87) 
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Here the function Kp(m 7r7T ) is defined by 

^p(m w ) = f 2 7r7r g P (m 7r7r ) + -(ln^ + - . (88) 

One sees that hard exclusive production of nir pairs at invariant masses close to the threshold = Im-^ opens 
a possibility to measure the coefficients L\\-\z of the EChL. These coefficients cannot be accessed experimentally in 
low energy processes as for such experiments one would need a strong source of gravity. Let us note that the couplings 
in the EChL involving gravity are also of relevance for the description of a meson gas out of thermal equilibrium ]39| ] 
(as e.g. in heavy ion collisions) as well as for the decay of a hypothetical light Higgs meson This simple example 
shows that hard exclusive processes have a big potential for studies of new properties of the EChL. In general, this 
kind of hard reactions allows to create low energy probes with quantum numbers that are not provided by nature or 
are hardly accessible. 

VI. RESULTS AND DISCUSSION 

In section || we have deduced analytical expressions for the cross sections and intensity densities in terms of skewed 
parton distributions of the nucleon, the electromagnetic pion form factor, the so called Omnes functions foim^) and 
/ 2 (m W7r ), and the momentum fraction M 2 carried by the quark in the pion. 

The results of section |l| show that M 2 enters the expressions for cross sections and intensity densities only in 
the combination 40M 2 Q /iV / + 2QM§/C F = 20/C F + (40/N f - 20/C F )M 2 Q . For N f = 3 this combination is only 
weakly dependent on the value of ■ Therefore we choose for the following analysis the asymptotic expression 
M 2 Q = Nf/(N f + AC F ) and state that our predictions are insensitive to the exact value of this parameter reducing 
in this way the number of parameters of our result^. Let us note that for this choice the contribution of the gluon 
27rDA to the isoscalar amplitude is exactly twice as large as the contribution from the quark 27rDA. 

For the numerical results shown in this section we use the models for the SPD's of the previous section. The pion 
form factor F^(m„) entering the isovector 27rDA is modeled using its parameterization given in p^| , which is in very 
good agreement with experimental data p3[ . The Omnes functions contain information about the nn resonances as 
well as about the non-resonant background. The function faim^") is dominated by the /2(1270) resonance resulting 
in a peak at m T7r = 1.275 GeV. In our analysis we use the expression ( |6l| ) for ^{jn^-^) with the D-wave pion scattering 
phase shift 6\ given by the Brcit-Wigner contribution of the /a(1270) resonance. For the the Omnes function /o(m Tff ) 
we use (except in Fig. [?]) the Pade approximation suggested in [Q . 

As a first result we show in Fig. the ratio of the differential cross sections for the production of pion pairs with 
isospin I — and / = 1 as a function of the invariant mass of the two pions m„ at the three different XBj -values 
^Bj = 0.1, XBj = 0.2, and XBj = 0.3. The squared momentum transfer t has been integrated from —0.6 GeV 2 to 
— *min = ~ m N x Bj/(^ ~ x Bj)- The plot shows that relatively large I — cross sections are to be expected at small 
TOtttt close to the threshold due to the S-wave contribution and around 1.3 GeV due to D- waves, which are dominated 
by the /2(1270) resonance. Near the p resonance around 0.77 GeV, however, the isoscalar contribution is negligible 
compared to the production of isovector pion pairs. The slight difference to our result in jl3| is due to the inclusion of 
the gluon SPD, which contributes only to isovector pion pair production, as well as to the fact that we have integrated 
over a finite Mnterval. 

Figure || shows that the most favorable kinematic region to observe the isoscalar channel by the measurement of 
intensity densities in 7r + 7r~ production is in the two m Tff regions above or below the p resonance and at relatively 
large values of the Bjorken variable xbj- 

This can also be seen from Figs. [| and |5|, where we show our results for the intensity densities (Pi(cos0)) 7r+7r and 
(P3(cos0)) 7r+Tr on the basis of ^-integrated cross sections as functions of XBj and m^^. The integration interval for 
the variable t is again as in the following figures —0.6 GeV 2 < t < — i m i n . Figure || shows that (P3(cos #)) 7r+7r is 
sizable only in the /2(1270) resonance region. This resu lts from the fact that this intensity density is proportional to 



the Omnes function f2(fn 7rw ) as can be seen in Eq. (A8) in the appendix 



4 Note that the process 7*7 — > tttv, which is also sensitive to gluon 27rDA's (see Ref. 0|), depends strongly on the ratio 
/M° as noted recently in Ref. Q 
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In Figs. ^|and|^ we show our results for the combination of intensity densities (Pi(cos0)) 7r+7r +y / 7/3(P3(cos6')) 7r+7r 
for t- and xbj- or t- and TO^-integrated cross sections respectively as functions of rn nn or XBj- The contribution of 
(P3(cos6')) 7r+7r in this combination is relatively small, especially in the region of small m^-values near the threshold. 
The choice of the kinematic region is motivated by the above results and by the kinematic range of the HERMES 
experiment, where corresponding measurements can be done f45f| . 

Let us mention that the expressions for the cross sections entering our results for the intensity densities depend 
on the variable y = p ■ q/p ■ I only through the common factor (1 — y)/y 3 in Eq. ( [f0| ) as long as we neglect the scale 
dependence of the SPD's. Therefore for comparison with experiments the cross sections can be taken as integrated 
over an arbitrary y-interval, a realistic choice could be 0.35 < y < 0.6. 

Especially in the region near the threshold m n7r — 0.28 GeV the predictions are very sensitive to the Omnes function 
/^(m^Tr). To illustrate this dependence we also show in Fig. ^the result for the alternative parameterization of p6] | for 
/ (m TW ) (dotted line). We see that the process of hard exclusive pion pair production is sensitive to the mechanisms 
of the low-energy scattering of pions in the isoscalar channel and therefore can be used to obtain new information on 
the chiral dynamics in this channel. 

We have checked that ignoring the D-term contribution to the quark SPD changes the results by about 10%. This 
may be taken as an estimate for the uncertainty of our predictions due to a lack of knowledge about the SPD's, 
especially about the functions E^r, £,i). Let us also note that NLO corrections in hard exclusive reactions, which are 
not included in the analysis of the present article, can be noticeable pljjlq] . 

So far we have restricted ourselves to electron-proton scattering. The corresponding results for a neutron target 
can easily be obtained by the interchange I~ <-> 1^ and Pj~ <-» 7j~ in Eqs. (^) and (|2^) due to isospin symmetry. A 
detailed analysis shows that in this case the ratio of the amplitudes for isoscalar and isovector pion pair production 
is approximately one order of magnitude smaller than for a proton target. This is due to a cancelation of the 
contributions of u- and d-quark 

SPD's in Eq. ( p^ ) because the u-quark density in the neutron (corresponding to the d-quark density in the pro- 
ton) is roughly half as large as the d-quark density. Consequently we can predict that also the intensity densities 
(Pi (cos 8))" *~ and (P 3 (cos0)) ,r 7r " 

are about one order of magnitude smaller for neutron targets. This observation 
might help to check the whole picture in experiments with deuteron targets. Such experiments could also help to 
disentangle the contributions of the different quark and gluon SPD's. 

1=0 , 1=1 
a /a 

0.4 
0.35 



0.3 
0.25 

0.2 
0.15 

0.1 
0.05 

0.0 

0.4 0.6 0.8 1.0 1.2 1.4 1.6 
m^GeV] 

FIG. 3. The ratio of the differential cross sections for isoscalar and isovector pion pair production at three different values 
for x-Bj as a function of m W7r . 
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FIG. 4. {Pi(cosO)) n * as a function of XBj and m W7T . 
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(P^cos 6») + (7/3) 1/2 P 3 (cos $)) 



+ - 

7T 7T 




0.6 0.8 1.0 1.2 
m^GeV] 



as a function of m n -n with cross sections integrated over xbj from 0.05 to 0.4. 



(P^cos fl) + (7/3) 1/2 P 3 (cos 0)) 



+ - 

7T 7T 



-0.02 



-0.04 



-0.06 




-0.08 



-0.12 



-0.14 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 



FIG. 7. (Pi(cos0)) ,r+lr + \/g(P3(cos 0)) 7T+7T as a function of ib j with cross sections integrated over m^^ from the threshold 
to 0.6 GeV. The dotted line shows the corresponding result obtained using the fit of instead of the Pade approximation for 
the Omnes function /o(m,r-/r). 
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VII. SUMMARY 



We have analyzed the angular distributions for exclusive electroproduction of pion pairs. Our analysis was based 
on the amplitudes at leading order in 1/Q 2 and a s . These amplitudes were expressed in terms of skewed parton 
distributions and two-pion distribution amplitudes, for which we have used realistic models. We have introduced 
intensity densities as Legendre moments of the two-pion angular distributions and we have shown that they are a 
good probe for the contribution of isoscalar pion pairs to the considered process. Our analytic results show that 
isoscalar pion pairs are mostly produced by two collinear gluons, which, in principle, opens a possibility to study 
the gluon content of the isoscalar irir states. We predict sizable effects that could be measured at experiments like 
HERMES. As a promising kinematic range for a corresponding experiments we identified the m^-regions near the 
threshold and around the /2(1270) resonance and relatively large values of the Bjorken variable XBj > 0.2. The m n7r - 
dependence should show a highly characteristic pattern. Our results are nearly parameter free such that experimental 
data would provide an excellent check for the whole formalism describing exclusive electroproduction processes with 
skewed parton distributions. In addition we have shown that studies of hard exclusive pion pair production near to 
the threshold open a new way to probe the effective chiral Lagrangian. 



ACKNOWLEDGMENTS 



We are grateful to M. Amarian, H. Avakian, B. Clerbaux, 
U.G. Meifiner, D. Miiller, P.V. Pobylitsa, S. Schaefer, E. Stein 
thanks J. Gasser for providing the details of Ref. 26 



L. Frankfurt, A. Kirchner, N. Kivel, L. Mankiewicz, 
j M. Strikman, and O. Teryaev for useful discussions. 
The work has been supported by the Russian- German 



M.V.P. 

exchange program, BMBF, DFG, COSY-Julich, and Studienstiftung des deutschen Volkes 



APPENDIX: 



In this appendix are shown the full expressions for the cos ^-integrated absolute squares of the T-matrix elements, 
weighted with Legendre polynomials Pi (cos 6), as they enter the intensity densities. 
Using Eqs. @ - M), @, and @ we get 



(Al) 
(A2) 



with 



and 



f dcos^|T^ +7r T= f dcos^(|T / =°| 2 + |T / = 1 | 2 ) 
J-i s , J-i s , 

f dcos^|T^°| 2 = f dco S 0V|T I=o | 2 

J-l g, J-l g, 

C dcosfl V it'=T - (^) 2 c| i T \ r 2 + 20*£V 

/-! V ' 18iV 2 p .l XBi yU 2l » '* { Nf + C F J 

{3 ° ~f f \Mm™)\ 2 + ^\h(m^)\ 2 



The other non- vanishing integrals are: 



pi 

/ doxe^T*** | 2 Pi(cosf?) = 

J — 1 ci/ 



(e7ra s ) 2 C 2 F 1 /40M 2 Q 20M 2 C 



18 A 2 p-lx^y\ N f C F 



(A3) 



£ dcosflg \T^? = p . ^ + ^ + i~/o\W\FAm^ ■ (A4) 



x2Re< 



2Al + IJ + -^Ig 
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3C-/3 2 , . „ 4/3 2 



J dcos6j2\T^ | 2 P 2 (cos0) = 



(e^a s ) 2 C| 



18 A^ 2 p-isnjl/ 



2/40M 2 Q 20M 2 G \ 2 



TV/ C* F 
3C-/3 2 2 „ „ tl O , /? 4 



90 



144/3 2 



f 1 

J — l ci 



C F I 35 

(e TOs ) 2 C 2 1 V | _ V40M 2 Q , 20M G ^ 2 
18iV 2 p./^BjyV' " d 



3C - /3 2 fl2 



90 



f 1 /? 4 

/3 Re|/ (r7l 7 r7r)*/2(TO7r 7 r)| + ^\M m ^n)\ 



f decs. V |T^-| 2 P 3 (cos.) = Jj^m^_(^t + 



/ dcos6»^ IT^ 71 | 2 F 4 (cose»)= / dcosflV |T 7r ° 7rO | 2 P4( c ° s6 

V-l e, J-l a , 



187V C 2 p./xBjy^ 



V40M 2 Q 20M 2 G \ 2 /3 4 . , |2 



\ JV/ C F J 315 



(A5) 



(A6) 



(A7) 



J^-/2(m 7r7r )*P 7r (m 7r7r ) S (A8) 



(A9) 
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